Scaling relation for determining the critical threshold for continuum percolation of 

overlapping discs of two sizes 
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We study continuum percolation of overlapping circular discs of two sizes. We propose a phe- 
nomenological scaling equation for the increase in the effective size of the larger discs due to the 
presence of the smaller discs. The critical percolation threshold as a function of the ratio of sizes 
of discs, for different values of the relative areal densities of two discs, can be described in terms 
of a scaling function of only one variable. The recent accurate Monte Carlo estimates of critical 
threshold by Quintanilla and Ziff [Phys. Rev. E, 76 051115 (2007)] are in very good agreement with 
the proposed scaling relation. 
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In recent years, there has been a lot of interest in study- 
ing continuum percolation, owing to its many applica- 
tions. For a review on continuum percolation, see 
Continuum percolation of overlapping objects of various 
sizes and shapes, spheres and discs 0, H, 0, H, 0, ellip- 
soids Q, plates [1], sticks @, oriented cubes ff(| etc., 
has been studied. In applications like the modeling of 
porous media, one of the most important parameters is 
the distance from percolation threshold, and several ap- 
proximation schemes have been proposed to determine 
the percolation threshold for different types of disorder. 

In this paper, we discuss the case of continuum per- 
colation of overlapping discs of two sizes in a plane. We 
propose a phenomenological equation for the increase in 
the effective size of the larger discs in the presence of the 
smaller discs. We check our theory against data on crit- 
ical thresholds by Quintanilla and Ziff The agree- 
ment is found to be very good. 

We consider a percolation model of a mixture of circu- 
lar discs of two sizes randomly placed in a plane. Con- 
sider a finite area S and randomly drop discs in S. The 
probability that a given small areal element dA contains 
the center of a dropped disc is ndA, independent of other 
discs. Once a center of the disc is chosen, it is assigned 
a radius R\ with probability /, and i?2 with probability 
(1 — /). We denote the ratio of radii R1/R2 by A. The 
number density of discs with radius R\ is then n\ = fn, 
and that of radius R% is n% = (1 — f)n. The total num- 
ber density of discs, irrespective of radius, is n = n\ +112- 
We propose an approximate formula for the critical per- 
colation threshold in terms of A and /. We express this 
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function of two variables in terms of the function £(A) 
which gives the correlation length £ as a function of the 
areal density A of single-sized discs. 

The earliest proposal for determining the critical 
threshold for overlapping discs was by Scher and Zallcn 
0. They noted that the total covered fractional area at 
critical threshold was nearly constant for a mixture of 
discs of different sizs, if the polydispersity of the mixture 
was small. However, if the polydispersity is large, and 
one takes discs with several different radii, the total cov- 
ered fraction at critical threshold can be made as close 
to one as we wish The ori gina l heuristic arguments 
have been made rigorous later [13 ] . 

We start by summarizing the qualitative arguments 
of [||. Let us assume, without any loss of generality, 
that the R\ < i?2- We consider the plane on which the 
smaller discs of radius R\ each have been thrown in ran- 
domly with n 1 discs per unit area. The areal density 
of these discs is then A\ = irRfni. Note that A\ is a 
dimensionless number giving the ratio of total area of 
discs thrown in to the area of the plane. In the case 
of percolation of discs of equal radii, the areal density 
of the discs at the percolation threshold is independent 
of the size of the discs. Let this critical value of A be 
denoted by A*. We assume that Ay is below critical 
threshold A*, and the small discs by themselves do not 
percolate. From numerical simulations, the value of A* 
is known quite accurately A* as 1.128085. The corre- 
sponding value of the covered area fraction is given by 
(f>* = 1 - cxp(-A*) « 0.6763475(5) [Il|. 

The two point correlation function, G(r), is defined as 
as the probability that two points at a distance r from 
each other, chosen at random, belong to the same cluster 
when only the smaller discs have been dropped. Below 
criticality, this decays exponentially with distance, i.e., 
G(r) ~ exp(— r/£i). And using simple scaling invariancc 
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of the problem R± — > ai?i , we have 

= fliff(Ai) 



(1) 



where the function g(x) determines how the correlation 
length varies with areal density, and is independent of 
R x . 

Now we throw in a single disc of the larger radius R 2 , 
and look at the cluster of discs that are connected to 
this single large disc. Then, each such cluster looks like 
a somewhat bigger fuzzy disc of size R 2 + AR 2 . Let us 
assume that the variation between different clusters may 
be neglected. This approximation is quite good if R2 3> 
£1, but less valid if R 2 /£,i is not so large. The percolation 
problem can then be considered as a percolation of these 
larger effective discs. The number density n 2 of these 
effective larger equal-sized discs of radius R2 + AR 2 that 
have to be dropped to reach criticality is given by 



n* 2 Tr(R 2 + AR 2 ) 2 = A* 



(2) 



We will consider this equation as the definition of AR 2 . 
In the simple approximation 



Ai? 2 « c a 



(3) 



was used, where c is some constant of order one. This 
gives the correct limiting behavior that for any initial 
density A\ of the smaller discs, the critical value of the 
areal density of larger discs ^2(^2) tends to A* as R 2 
tends to infinity, keeping A\ fixed. Also, the other limit 
when we keep A 2 fixed at any value below A* , and slowly 
increase A\ till we reach critical percolation, then the 
critical value of A*(Ri) to reach criticality tends to A* 
as R\ tends to zero (l2j . 

However, Eq.© strongly underestimates the value of 
AR 2 . Consider two discs of radius R 2 thrown in a sea 
of randomly dropped smaller discs of areal density A\. 
Call these discs 1 and 2 and, let the minimum distance 
between these discs be denoted by D (Fig. 1). We denote 
by Prob£,(l ~-» 2) the probability that there is a path of 
overlapping smaller discs between the larger discs, and 
they belong to the same cluster. Thus, Pro6u(l 2) is 
a measure of the connectivity correlations in the problem 
of percolation of single-sized discs. 

Clearly, Prob[>(l 2) is a decreasing function of the 
separation D, which will decrease exponentially from 1 
to 0, as D varies from to infinity. For large D, this 
decreases as exp(— £)/£i). The dependence of this on R 2 
comes from the fact that the prefactor of the exponen- 
tial would depend on R 2 . Also, for D comparable to £1, 
the -D-dependence can not be approximated well by a 
simple exponential. However, we can define an effective 
size Ai?2 ff by the requirement that this probability is a 
fixed value, say 1/2, when D = 2AR 2 . Then, a better 
estimate of Ai?2 than Eq. © is given by 



Ai? 2 » ARf . (4) 

The ARf as defined is a function of R\ , R 2 and ni 
(or, cquivalcntly £1). For D comparable to £1, we cannot 




FIG. 1: Two large discs of radius R2 in a background of ran- 
domly dropped smaller discs. The least separation between 
the discs is D. 



use the large D exponential decay of Pro&£>(l 2) to es- 
timate Ai?2 ff - However, if R\ <C £1, then we can assume 
that the leading dependence is from £1, and correction 
terms involving powers of can be neglected. Then, 

Ai?2 ff , to leading order, is only a function of R 2 and £x- 
Using the fact that the probabilities are invariant if all 
distances are scaled by same factor, we get 



ARf 



(5) 



where h(x) is some, as yet unspecified, scaling function 
of its argument x. Now, clearly, Probo(l 2) is a 
monotonically increasing function of R 2 , which tends to 
1 as R 2 tends to infinity, keeping D fixed, as then the 
problem is that of percolation in a very long strip, and 
somewhere or other, there will be a connection of smaller 
discs. This implies that ARf must tend to infinity if 
R 2 tends to infinity. Also, in the case fliC^Cdi it 
must tend to infinity as £1 tends to infinity. The simplest 
form of h(x) that is consistent with these requirements 
is a simple power-law form, which gives 



Am 



k i?2~ 



(G) 



Here k is some constant of order 1. The main im- 
provement in this form over Eq.Q is the inclusion of 
dependence on R 2 . 

The power-law dependence of R 2 a on R 2 is seen 
most easily by considering a perturbation expansion of 
Probo{^ 2) in powers of n\. Let 



Prob D {\ -w 2) = mFi(D, R 2 ) + n\F 2 (D, R 2 ) 



(7) 



In the first order in m, the configurations that con- 
tribute to Probn(l 2) are those where a single small 
disc overlaps with both the bigger discs. This is possible 
only if D < 2R\, and in that case, if there is at least 
one small disc in the region which is within a distance 
R\ + R 2 from the centers of the discs 1 and 2 ( See Fig. 
[2]). For small m, the probability of this event is propor- 
tional to the area of the shaded region in Fig. [5] Using 
elementary geometry, it is easily seen that for R 2 3> R\, 
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FIG. 2: Two large discs of radius R2 with separation between 
the discs D. A larger circle of radius Ri + R2 is drawn sur- 
rounding each of the discs. If the center of any small disc falls 
in the intersection region (shown shaded) of the larger circles, 
it forms a connecting path betwen them. 

the area is proportional to R^ 2 {2R\ — D) 3 ' 2 . Thus we 
get 

F X {D,R 2 ) ~ i?2 /2 (2i?!-L>) 3 / 2 , for < D < 2R X . (8) 

Thus, we see that F\, and by extension Pro6o(l 2) 
has a strong dependence on R 2 . Of course, for £1 ^> Ri, 
higher order terms in ri\ make significant contribution, 
and they would change the precise form of the functional 
dependence on R 2 . 

Again, we assume that the larger discs act as discs of 
radius R 2 + AR 2 , with AR 2 ~ ARf , given by Eq.©. 
Expressing n 2 in terms of A 2 , the areal density of the 
larger discs, the criticality condition may be written as 

—l = ^/A 2 -l^k[Xg(A 1 )] a (9) 

The above equation is clearly invariant under scaling 
of all lengths by the same factor. We can determine the 
value of a, in the limit £1 3> R 2 . Then, assume A\ ~ 
A* {I — e). Then, e < 1 implies that £1 > R±. 

Clearly, the number density of additional discs of ra- 
dius R 2 required to reach criticality would be less than 
with discs of size R\. Hence, in terms of areal densities, 
this bound becomes A*. < A*e\~ 2 . Also, as discussed in 
@ , the total areal density of discs at criticality is greater 
than A* when all discs are not of same size, A\ > eA*, 
Thus, A* 2 ~ e. Then, AR 2 ~ e -1 / 2 . Since it is known 
that g(x) ~ (A* - x)~ v for x near A*, with v = 4/3. 
Thus, comparing powers of e we see that a = 3/8. 

Our proposed approximation can be directly checked 
against numerical data. Quintanilla and Ziff have given 
a very extensive table of data giving different values of 



Ai,A 2 for different values of Ri/R 2 , that define critical 
surface [H. Using Eq.®, if we plot Y = \~ a [^ A* /A 2 — 
1] versus X = A* — A\, all points should fall on a sin- 
gle curve Y = g(A* — X) a . The result is shown in Fig. 
3, where we have plotted data corresponding to five dif- 
ferent values of A = 0.10,0.20,0.30 and 0.50 We get a 
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FIG. 3: Scaling collapse of the Monte Carlo data of [Tl|] . Re- 
scaled data of AR2 is plotted vs A* — Ai, the deficit in the 
areal density Ai of smaller discs from the critical value for 
mono-disperse discs A*, for different values of the ratios of 
radii A. 



very good collapse. We do not show other values, in or- 
der not to clutter up the figure, but have checked that 
the collapse is as good with them as well. Note that no 
free parameters have been used to generate the scaling 
collapse. 

Define $(21) = [<7(a;)] a . The function &(x), which gives 
the equation of the curve is in principle calculable if we 
can solve the problem of percolation probability with sin- 
gle sized discs. As of now, we only know the behavior of 
$ in certain regimes. For small x, ~ x and for x 

near A*, varies as (1 - x/A*)' 1 / 2 . Hence we pa- 

rameterize the curve as 

<f>(x) ~kx(l + cx){\- x/A*)- 1 ' 2 . (10) 

The values k = 0.25 and c = 2.20, give a fairly good fit. 
The curve using these fitting parameters is also shown in 
Fig. 2. 
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